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Abstract: 

Dispersion-managed solitons are reviewed within a Gaussian variational approximation and an 
integral evolution model. In the normal regime of the dispersion map (when the averaged path 
dispersion is negative), there are two solitons of different pulse duration and energy at a fixed 
propagation constant. We show that the short soliton with a larger energy is linearly (exponentially) 
unstable. The other (long) soliton with a smaller energy is linearly stable but hits a resonance with 
excitations of the dispersion map. The results are compared with the results from the recent 
publications ^. 
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1 Introduction 



New ways in optimization of existing telecommunication systems based on dispersion management 
technology attracted recently wide research interest from soliton-based groups (see reviews fl], ) • 
The main idea was to combine a high local group- velocity dispersion with a low path-average 
dispersion. The former feature results in the reduction of the four-wave mixing while the latter 
one reduces the Gordon-Haus timing jitter effects. When the path-average dispersion is small 
and normal, i.e. the defocussing segment in the fiber is dominant over the focussing one, a new 
phenomenon of branching of soliton solutions was discovered ||, |, |, |, 0. The sohton propagation 
in this regime is not supported by an uniform-dispersion optical fiber and seems to be one of the 
remarkable achievement of the dispersion management with sufficiently high local dispersion. 

The stability of branching soliton solutions in the normal dispersion regime was a subject of 
intense studies which lead to contradictory conclusions. Grigoryan and Menyuk announced the 
linear and nonlinear stability of both the branches while Berntson et al. conjectured instability 
of one of the branches |^]. 

In this paper, we intend to shed light on the complicated issue of existence and stability of soliton 
signals in the normal regime of the dispersion map. We find, in the small-amplitude approximation, 
that there exist two branches of soliton solutions for different levels of energy and different pulse 
durations at a fixed propagation constant. The short pulses with larger energy are proved to 
be linearly unstable, while the other (long) pulses with smaller energy are neutrally stable. We 
show that the transition from large-energy unstable solitons to the stable soliton signals occurs 
via long-term transient oscillations. The two branches of soliton solutions correspond to a single 
(small-energy) branch B in Fig. 1 of |p. Depending on a normalization condition (see Section 4), 
this branch may be either stable or unstable. 

Our strategy to develop the small-amplitude approximation is based on the combination of two 
analytical approaches: the Gaussian variational approximation and the integral evolution model. 

The Gaussian variational approximation, being inaccurate in details, is still useful for a quick 
and rough analysis (see Q for review and references). Also it was shown that the method can be 
extended to a rigorous Gauss-Hermite expansion of the basic model We improve the previous 
results summarized in by deriving a new dynamical system from the variational equations of 
a Gaussian pulse. The system clearly displays the linear and nonlinear instability of the short 
Gaussian pulse with larger energy. 

More rigorous analysis of the problem is based on the integral evolution model obtained by 
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Gabitov and Turitsyn |jT^ and by Ablowitz and Biondini |11|. Although this model is more compli- 
cated from computational point of view (see recent papers |12, |l3|), we managed to study numer- 
ically the construction of the linear spectrum of dispersion-managed solitons. Our results confirm 
the instability and transition scenarios predicted within the variational model. We also deduce from 
this model that the soliton signals in the normal regime of the dispersion map are in resonance with 
the wave continuum of linear excitations of the map. The resonance implies isually the generation 
of wave packets from stable pulses oscillating in time. The latter effects are beyond the accuracy 
of the analytical model and are left for further studies. 

2 Gaussian approximation: New dynamical model 

We study the NLS model in the dimensionless form Q, 

iuz + ^D{z)utt + e (^^DqUu + \u\'^u^ = 0, (1) 

where u{z, i) is the envelope of an optical pulse in the retarded reference frame of the fiber. The 
small parameter e measures the length of the dispersion's map and the inverse variance of the local 
dispersion. After normalization, Dq and D{z) are assumed to be of order of 0(1), and 

(D) = r D{z)dz = 0, D{z + l) = D{z). (2) 
Jo 

Further physical motivations for derivation and verification of the model (||) can be found in ^ . 
Soliton-like optical pulses are solutions of the model in the form, 

n(z,^)=V(^,^)e-'^^ (3) 

where n is the propagation constant and ip{z, t) is a soliton profile satisfying the boundary condi- 
tions, 

tP{z + l,t) = i^{z,t) (4) 

and 

Urn ^p{z,t) = 0. (5) 

One of the conventional approximation for soliton solutions of NLS-type equations is based on 
averaging the Gaussian anzatz in the Lagrangian density and further varying the Lagrangian density 
with respect to parameters of the Gaussian pulse (see for review). The Gaussian approximation 
is the first term of the Gauss-Hermite expansions when solving the NLS equation (||) in the limit 
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e — s- . We apply the Gaussian anzatz in the form, 

Here the four parameters of the Gaussian pulse are: c{z) - the amplitude, ^(z) - the gauge pa- 
rameter, a{z) - the pulse duration, and P{z) - the chirp. It was found that the four equations for 
variations of the Lagrangian density can be decoupled into a system for a{z) and P{z) of the form, 

da 4eEa^/^P 



dz (a2 + 4/j2)3/2' 



(7) 



i = + 2(.2+ 4/^^)3/2 J- 

The phase factor (j)[z) is expressed in terms of a[z) and 

while the amplitude c{z) is given in terms of the input energy constant E as 



= ^ [\z r dt\u\\z,t)>0. (10) 



\/2a a/vFJo ^-oo 

The stationary pulse (|3|) - (|5[) corresponds to the periodic solutions of the system (^) and (^) in the 
form, 

a{z + l) = a{z), p{z + l) = p{z), 0(z + 1) = </.(z) + e/i. (11) 

For simplicity, we study the periodic solutions in the limit e ^ by using a two-step dispersion 
map with zero average, 

r,f \ / ^1' 0<z <L 

^^") = \ L<z<l ' (12) 

where 

m = DiL = -1)2(1 - L)>0. 
The asymptotic solution in the limit e ^ can be sought in the regular form, 

aiz) = as + 0(e), f3{z) = D{z')dz' + /?, + 0(e), 

JO 

where Os, Ps are constant. The periodic solutions appear when Ps = —'m/2 and Os is a root of the 
equation. 



Do = Eg?J'^ 



log' 



(m^-l- 0^)1/2 2m \ as 



(13) 
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In addition, the propagation constant fj, can be obtained as a function of E and according to 
the equation, 



1 



5 , (m+irm? +a2)V2 
+ — log 



(14) 



These equations have been already derived in the literature, see [P, ^ for ( [Tsl) and |^] for (14). 
However, the relations (^) and (^4|) were viewed typically under the normalization condition. 



m 



E 



2S 



a. 



1 

5' 



(15) 



where S is called the map strength. In this normalization, the expression (|13| ) gives a small- 
amplitude limit of the results of |5|, i.e. the slope E/Dq is a function of S. The existence of 
solitons was identified both for Dq > (when S < Sthr) and for Dq < Q (when S < Sthr), where 
Sthr ~ 3.32. 

In this paper, we develop a different frame to view the soliton solutions ( [l3| ) and ([l^). Guided 
by the stability analysis of solitons in generalized NLS equations |14], we fix the parameters of 
the model {DQ,m) and construct periodic solutions as a one-parameter family in terms of the 
propagation constant /i. As a result, the parameters a and E can be found from ([T^)-([T^) as 
as = OsifJ-) and E = E{^). These functions are shown in Fig. l(a,b) for Dq = 0.02 and m = 2 and 
in Fig. 2(a,b) for Dq = —0.02 and m = 2. Obviously, the branching occurs at Dq < (i.e. at the 
normal regime of the dispersion map), when the dispersion map is defocussing on average. The two 
solutions coexist for fj, > fj,thr{E>Q,m) and E > Ethr{Do,m). Both the branches I and II correspond 
to a single branch B in the small-amplitude approximation under the normalization condition ( |l5|) 

i- 

In order to describe non-stationary dynamics of the Gaussian pulse near the periodic solutions, 
we derive a dynamical model from Eqs. (0) and (P) by setting, 

a(z) =ao(C) + eai(^,C) + 0(e2) 



and 



Diz')dz' + (3o{0 + e(3iiz,C) + 0{e^). 



Here ^ = ez is the distance to measure the evolution of a Gaussian pulse over many map's periods. 
The coupled system (|7|) and (^) can be averaged over the map's period subject to the periodic 
conditions: ai{z + 1,C) = oi(-2, C) and (3i{z + 1,C) = Pi{zX)- Then, the non-stationary system 
reduces to the dynamical model for ao(C) and /3(C), 



dao 



Ea, 



5/2 



(a2+4/?2)l/2 (a2+ 4(^^ + ^)2)1/2 



(16) 



dC 



Ff3iao,fio) = Do 



Ea 



3/2 



4m 



4(m + /3o) 



4/3o 



2/3o + (ao + 4/32)1/2 



(17) 



^2(m + po) + (ag + 4(m + PoW^^ , 
This system has of course the same stationary solutions ao = cts and /?o = /3s = —m/2 as those 
given in ([131) . The stationary solutions appear as equilibrium states in the dynamical system, whose 
stability can be found by linearizing, 

ao(C) = as + Aae'^^, 
/3o(C) = /3s + A/3e*^f, 



where the eigenvalue A is 



2Ea 



3/2 



3Do + 



at 



{m? + 0^)3/2 



(18) 



dp ^"""^ da ^'""^ (m2 + 0^)3/2 

We plot in Fig. 3 to confirm that > for branch I of the periodic solutions and A^ < for 

branch II (cf. Fig. 2). Thus, the linear analysis predicts the instability of the short Gaussian pulses 
with larger energy at a fixed propagation constant /i (branch II). In the limit /x — > fJ,thr{Do, rn), the 
instability disappears, i.e. 

lim X^ifi) = 0. 

This property follows from Eq. ([TtI) in the limit fi ^tfir{DQ, m), when 

da 



-dE^'^'^^'^dH^^ 



dfi 



0. 



(19) 



Connecting Eqs. ( [l^ ) and (p!9|), we find the following asymptotic approximation, 

^2 ^ r dFp/dE.dFjdp \ dE^ ^20) 
V da/d^ ) djjL 

Taking into account that dE/dfi, dFp/dE, and dFa/dp are all positive for Dq < 0, and a ~ (/u — 
P'thrY^'^ (see Fig. 2(a,b)), the asymptotic approximation (|20| ) produces the result, A ~ {^jl — nthrY^'^ ■ 

The nonlinear dynamics of the system ([T^ ) and ( |T7|) is shown in Fig. 4 for Dq = —0.02 and 
m = 2. At a fixed value of the energy E, there are two stationary Gaussian pulses of different 
durations: a short pulse is a saddle point, while a long one is a center. Inside the separatrix loop, 
there are oscillations of the pulse trapped by the center point. Outside the separatrix, the Gaussian 
pulse transfers to chirped linear waves. 

We notice that the transition scenario resembles the nonlinear dynamics of unstable solitons 
in generalized NLS equations [O]. The only difference is that the unstable branch in generalized 
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Figure 1: Energy E (a) and pulse duration as (b) for the Gaussian pulse (0) at Dq = 0.02, m = 2. 
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Figure 2: Energy E (a) and pulse duration as (b) for the Gaussian pulse (^) at Dq = —0.02, m = 2. 
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Figure 3: Growth rate for the Gaussian pulse at Dq = —0.02, m = 2. The branch II with 
< corresponds to the unstable Gaussian pulse. 
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Figure 4: Phase plane (/?, a) for the nonhnear dynamics of a Gaussian pulse (^) at Dq = —0.02, 
m = 2 and fj, = 1. 



NLS equations is located for those values of soliton propagation constant /x, where dE/d^ < 0. 
Although this conventional stability criterion failed for the dispersion-managed solitons (see Fig. 
2(a)), the instabihty development shows up to be alike (cf. Fig. 4 here and Fig. 2(b) in p^]). 

3 Integral evolution model: Numerical analysis 

The Gaussian approximation of the optical pulse in the NLS model (|^) can be improved by sum- 
mating all higher-order Gauss-Hermite solutions of the linear equation, iu^ + 0.5D{z)utt = as 
shown in However, this analysis results in a complicated infinite-dimensional system of alge- 
braic equations for parameters and coefficients of the Gauss-Hermite expansion. Instead, we adopt 
a direct asymptotic method to deduce an integral evolution model valid in the limit e ^ 0. This 



method is based on Fourier expansion of solutions of the linear equation above [|^, 11] as well as 
on the asymptotic expansion, 

u{z, t) = uo{z, t) + eui{z, t) + O(e^), 

where U(){z,t) is given in the Fourier form as 

ito = ^ y"°° dujW{uj, C) exp (^-^oj'^ (yj^ D{z')dz'^ + iuj?j . (21) 

Here W{loX) is a complex Fourier coefficient and ( = ez is the distance to measure the pulse 
evolution over many map's periods. By supplying the periodic condition ui{z -|- l,t) = ui{z,t) in 
the Fourier form, the NLS equation (|l|) can be reduced to the integral evolution model. 



iW(^ - -DoUJ''W + I I duJiduJ2r{uJiuJ2)W{uj + uJi)W{uj + uj2)W{uj + uJi + UJ2) = 0, (22) 



where 



r(uiiUJ2) = [ dzexp (iujiuj2 [ D(z')dz' 
47r^ Jo V ^0 

For the two-step dispersion map ([l2|), the integral kernel r(x) can be computed explicitly as 

1 sin(^) 



It is obvious that the dynamical system (|16|) and (|17D studied in the previous section can be 
deduced from ( p2|) within the same Gaussian approximation. This correspondence implies that the 
qualitative results on instability of short Gaussian pulses for Dq < can be reconfirmed within a 
more systematic theory. 

In this section we present numerical results consisting of three subsections. In the first sub- 
section, we construct a numerical solution of the stationary problem identifying optical solitons 
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in the normal regime, when Dq < 0. In the second part, we analyze the linearized problem and 
locate numerically the linear spectrum in the problem, indicating possible instability of optical 
solitons. Then, we simulate numerically the non-stationary problem described by ( [2^ ) and discuss 
the transformation routes of the unstable dispersion-managed solitons. 

3.1 Stationary solutions 

The periodic-type localized solutions of the NLS equation in the form (0) are equivalent to stationary 



solutions of (22) in the form, 

W{io,C) = <^{io)e'''^, (24) 



where ^{uj) is the real Fourier coefficient which defines ip{z,t) according to ([2l|). This function 
satisfies a nonlinear integral boundary- value problem, 

Ii+]^Dqu?^^{uj) = R{uj) = J J duJiduJ2r{uJiuj2)^i^^ + + ^2)^i^ + ^^1 + ^^2), (25) 

where ^{—lo) = ^(uj) (the symmetry condition) and lim^;^oo ^(^^j) = (the boundary condition). 
For numerical analysis, we intend to use the Petviashvili's iteration scheme |T^: 

$W(^) ^ (^) 

for n = 0, 1, 2, .... Within this scheme, the right-hand-side R{uj) can be approximated at the n-th 
approximation by <I>'^"'-*(lij) provided a certain stabilizing factor is introduced for convergence (see 
(|27| ) and (|2^ ) below). However, the numerical scheme breaks down for Dq < due to resonances 
at LO = zizujresj where 

^res = ]j (26) 

Indeed for uj = ±uJres, the left-hand-side of (^) vanishes. [Here we notice that /i > for the 
Gaussian pulse solutions (^) of the NLS model (||).] In order to avoid resonances in the numerical 
scheme, we add and subtract a dummy positive dispersion term 0.5|Do|'^^*^'('-^) to the left-hand-side 
of (p5|). As a result, the scheme converts to the following map, 

where Sn is the Petviashvili's stabilizing factor given by 
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The factor Sn is unity at the stationary solution and serves therefore as an indicator for termination 
of the iterating procedure. We stop iterations when — 1| < 10^^. 

To use the map (p7[), we apply the Simpson's integration method, reducing complexity due to 
the symmetry: ^{—uj) = ^(uj). As a starting solution, the profile <I>(lj) can be approximated by 



the Gaussian pulse with parameters corresponding to the periodic solution (|T^) and (14) 



7r£^v2asexp ( —-agUj' 



(29) 
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Table I. Iterations of the stabilizing factor Sn for Dq = 0.02 and Dq = —0.02 (branches I and II). 



Table 1 shows iterations for the stabilizing factor Sn in the three different cases: (i) Dq = 0.02, 
(ii) Dq = —0.02 (branch I), and (iii) Dq = —0.02 (branch II). For all the cases, the other parameters 
are fJ- = 1 and m = 2. In the first case, the convergence is monotonic and the profile for stationary 
soliton ^(oj) is shown in Fig. 5(a). The numerical value for energy of the stationary soliton is 
shown in Fig. 1(a) by a bullet. In the second case, the iterations converge slowly to the stationary 
soliton shown in Fig. 5(b). Sometimes, the convergence is accompanied by a single oscillation of 
Sn near unity. The numerical value for the energy is shown in Fig. 2(a) by a bullet at branch 
I. In the last case, however, the iterations oscillate and finally diverge. Inspection of the profile 
$^"^(a;) at a final iteration shows that the iterations change the initial pulse drastically leading to its 
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(a) 



(b) 




Figure 5: The profile of the stationary pulse $(u;) in the boundary-value problem (^) at Dq = 0.02 
(a) and Do = —0.02 (branch I) (b). The other parameters are: m = 2 and fi = 1. The dotted line 
displays the Gaussian pulse (|29|) for the same parameter values. 



disappearance. Two conjectures follow from this fact. Either the shorter soliton with larger energy 
at branch II does not exist as a stationary solution of (^) or it is unstable within the iterational 



scheme (27). Since the short Gaussian pulse does exist (see Fig. 2), we incline to work along the 
second conjecture. The iterational scheme ( [27| ) is not relevant for the time-evolution problem and 
rigorous analysis of linearized problem is needed to confirm predictions of the instability of the 
short stationary pulse. 

3.2 Linear spectrum 

There are several forms of the linear problem associated to the NLS-type equations. We will use 
the matrix form which was studied in our previous paper [l^ subject to certain simplifications. 



(30) 



The matrix form appears upon perturbations of the stationary solutions of (22) as 

The vector w(a;) = {wi,W2)'^ can be shown to satisfy the matrix linear problem, 

Aw(tj) = - + -Dqw^^ (T3w(a;) + J duJi {2Ki{uj,uJi)a3 + K2{uj,uji)a2) '^{ui), 

where 



(31) 
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Figure 6: (a,b) The linear spectrum A of Eq. ( |3l| ) for the stationary pulse at Dq = 0.02, m = 2 
and n = 1 which corresponds to Fig. 5(a). 

and the integral kernels are 

/oo 
duj2r[{uj - uji){lj - u;2)]<J> (^2)^(1^1 + uj2 - w), 
-00 

/oo 
duJ2r[{ui2 - UJl){uJ2 - UJ)]^ {i02)^{u!l - UJ2+ Uj). 
-00 

Using the stationary solution ^(uj) from the previous subsection and implementing the Simpson's 
integration method again, we solve the linear problem by using the linear algebra packages built 
in Matlab 5.2. We identify two types of modes of the linear spectrum: symmetric eigenfunctions, 
when w{—uj) = w(i^), and anti-symmetric eigenfunctions, when w{—uj) = —w{uj). 

The linear spectrum for the soliton of Fig. 5(a) is shown in Fig. 6(a,b). It consists of three main 
parts: (i) continuous spectrum, (ii) neutral (zero) modes, and (iii) internal (oscillatory) modes. 

When Dq > (the anomalous regime of the dispersion map), the continuous spectrum is 
located at the real axis for |A| > /i (see Fig. 6(a,b) where ^ = 1). Indeed, the continuous modes 
are singular in the Fourier representation, i.e. w(ti;) ~ 5{uj — Q). Then, the linear problem ( pl[ ) has 
the continuous spectrum at X = ±Xq, where 

\n = fi + ^Don^, (32) 

provided the following integral kernels are not singular, 

1 

47r^ J-00 

/oo 
~oo 
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Figure 7: (a,b) The linear spectrum A of Eq. ( |3l] ) for the stationary pulse at Dq = —0.02, m = 2 
and n = 1 which corresponds to Fig. 5(b). 

The neutral (zero) modes always appear at A = as double degenerate states for both symmetric 
and asymmetric eigenf unctions. However, the inaccuracy of the numerical method destroys the 
degeneracy of the zero modes. As a result, the two zero modes may appear either for small real or 
for small imaginary values of A. Fig. 7(b) displays two imaginary eigenvalues of order of O(10~^) 
which appear to be shifted from the origin of A due to this numerical effect. Since the zero modes 
are not of interest from stability point of view, we neglect this effect and leave the scheme without 
any additional modification. 

The internal (oscillatory) modes are located in the gap of the continuous spectrum as < |A| < 
fi. The set of internal modes may contain different number of eigenvalues. We have shown in 
previous paper that the set is empty in the NLS limit (which corresponds to the limit /i ^ at 



fixed Dq > and m). [Note that in |16| we used the Gaussian pulse (29) for approximating ^{uj) 
while here we substitute the numerical result from Eq. (p5|).] Then, we showed that the number 
of internal modes increases with the map's strength (if the parameter as is set to unity, the map's 
strength is proportional to m and vice versa). In Fig. 6(a,b) for /u = 1, we identify 14 internal 
modes for symmetric eigenfunctions and 12 internal modes for anti-symmetric eigenf unctions. Still 
complex eigenvalues are absent for Dq > which confirms stability of dispersion-managed solitons 
in the anomalous regime. 

The linear spectrum for the soliton of Fig. 5(b) is shown in Fig. 7(a,b). The continuous 
spectrum is seen to have changed drastically. When Dq < (normal regime of the dispersion map), 
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Figure 8: (a,b) The linear spectrum A of Eq. (|3l|) for the stationary pulse at Dq = —0.02, m = 
and = 1 which corresponds to the Gaussian pulse at branch II. (c,d) The symmetric (c) an 
antisymmetric (d) eigenvectors w{uj) for the unstable eigenvalues with Im(A) > 0. 
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the continuous spectrum covers the segment |A| < /x twice according to (32). As a result, neutral 
and internal modes, if any, become embedded in the wave continuum as seen in Fig. 7(a,b). This 
indicates a resonance of stationary soliton with the linear spectrum in the normal regime of the 
dispersion map. However, this resonance does not result to any instability of solitons of branch I 
within the linear theory. We discuss the resonance issue in Section 4. The two imaginary eigenvalues 
on Fig. 7(b) appear from the origin as artifacts of the numerical scheme as it was explained above. 

At last, we would like to construct the linear spectrum for the solitons of branch II. However, 
the stationary solutions were not identified within the Petviashvili's numerical method. Therefore, 
assuming that the solutions still exist, the profile ^(k) can only be approximated by the Gaussian 
pulse ( p9|) as we did in jl^. The linear spectrum in this approximation is shown in Fig. 8(a,b) 
for Dq = —0.02, m = 2, and /u = 1. The same type of the continuous spectrum is clearly seen not 
to possess any gap in the origin. In addition to this, we identify new complex eigenvalues both 
for symmetric and anti-symmetric eigenf unctions. These complex eigenvalues has relatively large, 
order of 0(10^^), imaginary part and they are associated with the instability of the solitons of 
branch II. The numerical result for the instability eigenvalues is in a reasonable comparison with 
the asymptotic predictions which follow from the dynamical model (^) and ( [T7| ) (see Fig. 3, branch 
II). The eigenvectors w{uj) for the unstable (imaginary) eigenvalues are shown in Fig. 8(c,d) for 
the symmetric and anti-symmetric eigenfunctions, respectively. The numerical approximations of 
the eigenvectors, being inaccurate in details, display clearly that the unstable modes are localized 
at the intermediate wave frequencies lo of the pulse spectrum (at lo ~ 4.5). Thus, the development 
of the unstable eigenvectors would affect the duration of the soliton pulse in the nonlinear stage as 
described in the next subsection. 

3.3 Non-stationary evolution 

To confirm the transition scenario, we simulate the non-stationary dynamics of unstable solitons in 
the integral model (|2^ ) by using the central-difference scheme, 

HtxT — = M / d^ic^^2r(L^it^2)e'^»"^"^^"y("H^+^i)^^"H^+^2)^^"^(a;+c^i+u;2) 

(33) 



where 



yW(^) = ty(cu,C„)exp r-Dou^Cn 



and Cn = nA^. An initial iteration can be done within a forward scheme starting with the initial 
Gaussian pulse (pi) with the parameter and E corresponding to branches I and II at = 1 on 
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Fig. 2(a,b). Evolution of a stable long pulse (branch I) is shown in Fig. 9(a,b), while that of an 
unstable short pulse (branch II) is shown in Fig. 9(c,d). The stable long Gaussian pulse quickly 
transits to the stationary pulse given by Eq. ( [2^ ) (Fig. 5(b)) which propagates later without visible 
distortions. The second and third peeks in the signal spectrum and profile (Fig. 9(a,b)) appear 
in complete agreement with the profile of the stationary DM soliton (cf. Fig. 1 from [|11|). Some 
oscillations along the distance ( are excited due to the difference between the Gaussian pulse ( p9D 
and the exact stationary solution of Eq. ( p5|) . These oscillations are small compared to the soliton 
profile and they do not change the duration of the soliton pulse (see Fig. 9(a)). The nonlinear 
resonance at uires i^res = 10 for Fig. 9(a)) is not seen to be excited during the signal propagation. 

Evolution of the short (unstable) pulse differs drastically from the previous picture. The pulse 
is being broaden during the evolution, it generates the strong radiation tail and tends to the long 
(stable) signal which has the first node at w ~ 4.5 (cf. Fig. 9(a) and Fig. 9(c)). This transformation 
is accompanied by the intermediate oscillations around the soliton's shape. Thus, we confirm the 
analytical predictions that the short pulses are linearly unstable and switch into long stable solitons 
via long-term oscillatory dynamics. The unstable eigenvectors at Fig. 8(c,d) clearly match at Fig. 
9(c) with the growing deformations of the localization of the pulse spectrum. 

4 Discussion: Resonance of dispersion-managed solitons 

We have shown that both the Gaussian variational approximation and the integral model prescribe 
the instability of short nonlinear signals in the normal regime of the dispersion map. This instability 
broadens the signal's profile through some intermediate oscillations. The long signals propagate 
stably then. 

We point out that the small-amplitude approximation considered here corresponds to the asymp- 
totic limit E/Dq (see branch B at Fig. 1 in Q). The stability of the small-amplitude branch 
for Dq < was previously under discussion in the literature [^, |5|, |l^. This discussion can be 
resolved if one takes into account the normalization condition (|l5|) , which was used in the previous 



works. We have checked that the normalization condition (15) selects the pulse solution along the 
stable branch I for Sthr < S < Sstab^ where Sgtab ~ 6.76. For S > Sstab, it selects solutions along the 
unstable branch II (see Figs. 2(a,b)). Thus, for the intermediate map strengths (when S < Sstab), 
the soliton signal propagates stably in the limit of small energies, as reported in Q. However, if 
the map strength S exceeds the value Sgtab^ the soliton signal breaks down and switches to a longer 
signal, as conjectured in 0]. 
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(c) 



(d) 



Propagation Distance =19.8 



Propagation Distance =19.8 




FREQUENCY DOMAIN 



TIME DOMAIN 



Figure 9: Propagation of the stable (a,b) and unstable (c,d) Gaussian pulses in the model (p^) at 
Dq = —0.02 and m = 2. The dotted line displays the initial pulse (p9|) for the branch I (a,b) and 
for the branch II (c,d) at // = 1. The solid line displays the profile at C = 20. 
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For both the branches, we observe the resonance appearing between the hnear wave spectrum 
and the stationary signals. This resonance is related to the fact that the origin A = is absorbed in 
the continuous spectrum of the linear problem. Another way to find the resonance is to construct 
the linear spectrum for the integral model (22), W{u!,C) ~ Woe*^^'^^'', where 

n = -^Dou^ > 

and Dq < 0. Therefore, for any ^ > there exists uJres such that Q{iJres) = A*- This resonance does 
not show up in the linear theory since the discrete and continuous spectra are separated. However, 
in the nonlinear stage, the resonance generally leads to emission of wave packets and soliton's decay. 

Since the transformation or decay of long stable solitons for Dq < have never been observed 
numerically (nor in our simulations reported in Figs. 9), it is likely that the effective gap in the 
spectrum still appears in the nonlinear theory. Also the truncated approximation given by the 
integral model ( |2^ ) may not be valid for a correct description of the resonance. The latter issues 
remain open for further analytical consideration. 
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